Statistical coalescence model 
of J/ip production at the SPS and RHIC 



Mark I. Gorenstein 



Bogolyubov Institute for Theoretical Physics, Kiev, Ukraine 
Institut fur Theoretische Physik, Universitat Frankfurt, Germany 



Abstract 

A recently developed statistical coalescence model of the J/ if) produc- 
tion is presented. The NA50 data on the J/tjj suppression pattern and 
the transverse mass spectra in Pb+Pb collisions at the SPS are ana- 
lyzed. The model predictions for the RHIC energies are formulated. 
The measurements of J/ip in Au+Au collisions at RHIC are crucial for 
disentangling the different scenarios of charmonia formation. 

I Introduction 

The experimental studies of nucleus-nucleus (A+A) collisions at the SPS and 
RHIC provide a rich information on hadron observables (multiplicities and 
momentum spectra). An extensive experimental program is motivated by a 
possibility to create a new state of matter - the quark gluon plasma (QGP) - 
in a laboratory. Any theoretical description of the QGP requires a macroscopic 
approach: statistical mechanics to describe the static properties of the matter 
and hydrodynamics to describe the system evolution. The equilibrium hadron 
gas (HG) model describes remarkably well the hadron multiplicities measured 
in Pb+Pb and Au+Au collisions at top SPS (y/s = 17 GeV per nucleon pair) 

II and RHIC (y/s = 130 GeV) @ energies, where the creation of QGP is 
expected. 

The HG model assumes the following formula for the hadron thermal mul- 
tiplicities in the grand canonical ensemble (g.c.e.): 
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where V and T are the system volume and temperature, respectively, rrij and 
dj denote particle masses and degeneracy factors. The particle chemical po- 
tential fij in Eq.(|l|) is denned as //j = bjfi B + Sjfig + c j^c > where bj,Sj,Cj 



denote the baryonic number strangeness and charm of particle j. The bary- 
onic chemical potential //g regulates the baryonic density of the HG system 
whereas strange fis and charm \xc chemical potentials should be found from 
the requirement of zero value for the total strangeness and charm in the system 
(in our consideration we neglect small effects of a non-zero electrical chemical 
potential). The total multiplicities Nj ot in the HG model include the resonance 
decay contributions: 

N tot = N . + Y J Br{R-^ 3 )N R , (2) 

R 

where Br(R — > j) are the corresponding decay branching ratios. The hadron 
yield ratios Nj 0t /Nf ot in the g.c.e. are the functions of T and and are 
independent of the volume parameter V. 

The temperature parameters extracted from fitting the multiplicity data are 
approximately the same for the SPS and RHIC energies: T = 170 ± 10 MeV 
[|], |J. This value is close to an estimate of the temperature T c for the QGP- 
HG transition obtained in Lattice QCD simulations at zero baryonic density 
(see e.g. 0). One may therefore argue that QGP formed at the early stage 
of A+A reaction hadronizes into a locally equilibrated HG, and the chemical 
composition of this hadron gas is weakly affected by the subsequent hadron 
rescatterings. 

2 Statistical coalescence model of J/ip production 

The production of charmonium states J/ip and ip' have been measured in A+A 
collisions at CERN SPS over the last 15 years by the NA38 and NA50 Collab- 
orations. These studies were mainly motivated by the theoretical suggestion 
[f§] to use J j ip meson as a probe for deconfinement in A+A collisions. 

Recently the thermal model |j] and the statistical coalescence model (SCM) 
[!) ( see a l so B H) H3) f° r the charmonium production in A+A collisions were 
formulated. The total J/ip multiplicity, Nj°L, in the thermal model || is given 
by Eq. (||), where Nj/^, N^>, N Xl , N X2 are calculated according to Eq.([l]) and 
Br(ip') = 0.54, Br(xi) — 0.27, Br(x2) — 0.14 are the decay branching ratios 
of the excited charmonium states into J/ip. The thermal model predicts 
that at high collision energies the J/tp to n ratio is independent of and 
the number of nucleon participants N p . This is because both (J/ip) and (n) 
multiplicities are proportional to the system volume and the hadronization 
temperature is expected to be approximately constant at high collision ener- 
gies. 



The SCM H assumes that charmonium states are formed at the hadroniza- 
tion stage. This is similar to the thermal model ||. However, in the SCM the 
charmonium states are produced via a coalescence of c and c (anti) quarks cre- 
ated by the hard parton collisions at the early stage of A+A reaction. The 
number of created cc pairs, N C c, in hard parton collisions, differs in general 
from the result expected in the equilibrium HG. One needs then a charm 
enhancement factor 7 C ||[] to adjust the thermal HG results to the required 
average number, (N C c), of cc pairs. The open charm hadron yield is enhanced 
by a factor 7 C and charmonium yield by a factor 7^ in comparison with the 
equilibrium HG predictions. This leads to a difference between the thermal 
model and SCM predictions for charmonia multiplicities. 

The canonical ensemble (c.e.) formulation of the SCM is M: 

= \ ■<< N ° 7WT + ^ N » • (3) 

where Nh and Nq are the HG multiplicities of all particles with hidden and 
open charm, respectively. The canonical suppression factor I\/Iq in Eq.(H) is 



due to the exact charm conservation (see e.g., [12], Iq and I\ are the modified 
Bessel functions). 

If ( N C c) is known, Eq. (H) can be used to find the charm enhancement factor 
7 C and calculate then the J /if) multiplicity: 

(J/VO = 7 C 2 N% , (4) 

where Nj°L is the total (thermal plus excited charmonium decays) HG J/ if) 
multiplicity. 

Note that the second term in the right-hand side of Eq. (D) gives only a tiny 
correction to the first term, i.e. most of the created cc pairs are transformed 
into the open charm hadrons. Eqs. (|3|,|4]) lead then to: 

W) = ^yfp . W « 1 ; (5) 

jytot 

W) = (N*) 2 -^-^ , (iV c - c ) » 1 . (6) 

Eq.@ assumes an exact conservation of N c — N 5 = and 'statistical fluc- 
tuations' of N C c = (N c + N s )/2 numbers. A more accurate treatment with 

1 This is formally analogous to the introduction of the strangeness suppression factor 
7 S < 1 El in the HG model. 



'dynamical' (Poisson-like) distribution of leads to ||: 

j^tot 

(JM = (N c5 ) (1 + (N c5 )) . (7) 

Eq.(0) coincides with Eqs.([|D at both (A^) << 1 and (A" cS ) » 1 limits 
(see Eqs.([5]||)). This is because at (N C c) « 1 the probabilities to create zero, 
P(0) = 1 — (N C c), and one, -P(l) — (A^), cc pairs are only important. On the 
other hand, at (N^) » 1 the statistical and dynamical fluctuations of N C c 
both obey the Poisson law distribution. 



3 J/ip suppression in Pb+Pb collisions at the SPS 

The centrality dependence of (N c5 ) in A+B nucleus-nucleus collisions can be 
calculated in Glauber's approach, (N C c)(b) = a^ 5 N TA B (b), where b is the impact 
parameter, Tab (6) is the nuclear overlap function and a^ N is the cc produc- 
tion cross section for nucleon-nucleon collisions. As discussed in Ref. 0, the 
deconfined medium can substantially modify charm production at the SPS, 
i.e., a^ 5 N in A+B collisions can be different from the corresponding cross sec- 
tion measured in a nucleon-nucleon collision experiment. The present analysis 
treats a^ N at the SPS energy as a free parameter. Its value is fixed by fitting 
the NA50 data. 



In the NA50 experiments [|l]| the Drell-Yan muon pair multiplicity (either 
measured or calculated from the minimum bias data) is used as a reference 
for the J /ip "suppression pattern". The number of Drell-Yan pairs is also 
proportional to the number of primary nucleon-nucleon collisions: (DY')(b) = 
a^p,TAB(b), where ct^y' * s the nucleon-nucleon production cross section of 
fi + H~ Drell-Yan pairs. The prime means that the pairs should satisfy the 
kinematical conditions of the NA50 spectrometer. As the Drell-Yan cross 
section is isospin dependent, an average value is used: u^y' = a DY' I (AB). 
For the case of Pb+Pb collisions, A = B = 208 and a%f b = 1.49 ± 0.13 fjb 
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The quantity to be studied is the ratio 

n Dr(b) s (8) 
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B'flf = 0.0588 is the decay probability of J /if) into /i 4 'fjT . Only the fraction r] 
of /x + /i~ pairs satisfying the kinematical conditions of the NA50 spectrometer 
can be registered. We treat rj in Eq.(|8]) as one more model free parameter. 

In the NA50 experiment, the neutral transversal energy Et of produced 
particles was used to measure centrality of the collisions. This variable, how- 
ever, provides a reliable measure of the centrality only if it does not exceed a 
certain maximum value: E? < 100 GeV (see Ref. [§]). The average number of 
participants is a linear function of the transversal energy, N p = Ex/q, in the 
domain E T < 100 GeV (q 0.274 GeV). At E T > 100 GeV the number of 
N p does not change essentially as E T grows. Therefore, the data at E T > 100 
GeV do not represent a centrality dependence of the J / if) suppression pattern 
but rather its dependence on fluctuations of the stopping energy at approxi- 
mately fixed number of participants. The influence of such fluctuations on J / if) 
multiplicity can be studied in the framework of the SCM model WA . 
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Fig. 1: The dependence of the J/i/j to Drell-Yan ratio on the transversal energy. The normal 
nuclear suppression curve is obtained at <r a bs = 6.4 mb. The SCM lines are calculated using 
Eq.(|§|). The vertical line shows the applicability domain of the SCM, N p > 100. 



The SCM fit of the NA50 data with Eq.( 
for details). The model parameters are: 



is shown in Fig. 1 (see Refs. [§, [14 
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with T and /xg fixed by hadron yield data. The model requires rather large 
enhancement for the open charm production (up to a factor of about 3.5 within 
the rapidity window of the NA50 spectrometer). A direct measurement of 



the open charm in Pb+Pb at CERN SPS by NA60 would allow to test this 
prediction. 

The SCM does not describe the NA50 data at small values of Et (the 
peripheral collisions with E < 27 GeV). This can be also seen from the if>' 
data. For T = 170 MeV the value of the thermal ratio is (ip')/{J/ip) = 0.04. 
This is in agreement with data in Pb+Pb collisions at SPS for N p > 100, but 
is in contradiction with data in p+p, p+A and very peripheral A+A collisions. 
This fact was first noticed in Ref. |T5| . 



4 J /ip enhancement in A+A Collisions at the RHIC 

The number of directly produced cc pairs at the RHIC energies can be es- 
timated in the pQCD approach and used then as the input for the SCM. 
The pQCD calculations for cc production cross sections were first done in 
Ref. fl6| . For the cross section a(pp — > cc) of the charm production in p+p 



collisions we use the results presented in Ref. [17 1. This leads to the value of 
a(pp — > cc) = 0.35 mb at y/s = 200 GeV and the A/i-dependence of the cross 
section for yfs = 10-4- 200 GeV is parameterized as |i0| : 




a(pp - cc) = (7 • 1 - > ( 10 ) 



with a 3.392 /xb, M = 2.984 GeV, a 8.185 and (3 1.132. 

The number of produced cc pairs in A+A collisions is proportional to the 
number of primary N+N collisions, N^, which in turn is proportional to N^ 3 



IS 



W = N^(N p ) a{PP Z Ct) - C a(pp cc) < 3 , (11) 

<7 NN 

where cr™fi — 30 mb is the inelastic N+N cross sections, C = 11 barn -1 . 

The results of the SCM can be studied analytically according to Eqs.(|],|]) 
in the limiting cases of small and large numbers of (N&). For N C c << 1 one 
finds: 

" (JVcc) Nl V (!) ^ > ' • { 1 

where we use the energy dependence of the pion multiplicity per nucleon par- 
ticipant (ir)/N p oc (v^) 1 ^ 2 IHH which approximately works in the SPS-RHIC 



energy region. The behavior ([12]) corresponds to the J/ip suppression: the 
ratio R decreases with increasing of both y/s and N p . This takes place at 



the SPS: this energy is still too "low" as (N C c) < 1 even in the most central 
Pb+Pb collisions. However, the behavior of the J /if} to N C c ratio is changed 
dramatically at the RHIC energies (see Fig. 2 and Ref. [10 for details). 
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Fig. 2. Left: The energy dependence of the (J/tp) to (N c ^) ratio in central Au+Au collisions. 
Points are the predictions of the SCM for the RHIC energies: y/s = 56, 130, 200 GeV. 
Right: The iVp-dependence of the (J/tp) to (N C c) ratio. The lines are the predictions of the 
SCM. The dashed line corresponds to yfs = 56 GeV (the J/tp suppression), the solid line 
corresponds to = 200 GeV (the J/tp enhancement). 

In central Au+Au collisions at \fs = 200 GeV the expected value of (N C c) is 
essentially larger than unit. For (N C c) » 1 one finds (f3 = 1.1): 



R = 



(N C c) 



N, 
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(V~s)^ 1/2 iV p 1/3 • (13) 



5 The m^-spectra of J/ip and ip' and QGP hadronization 

In Ref. ||20|| we formulated the hypothesis that the kinetic freeze-out of J/ip and 
if)' mesons takes place directly at hadronization. The effect of rescattering in 
the hadronic phase was recently studied within a "hydro + cascade" approach 
[ f2~l| , ^] . A+ A collisions are considered there to proceed in three stages: hydro- 
dynamic QGP expansion ("hydro"), transition from QGP to HG ("switching") 
and the stage of hadronic rescatterings and resonance decays ("cascade"). The 
switching from hydro to cascade takes place at T = T c , where the spectrum 
of hadrons leaving the surface of the QGP-HG transition is taken as an input 
for the subsequent cascade calculations. The results of Refs. |2l|, ^] suggest 



that the transverse momentum spectra of Q is only weakly affected during the 
cascade stage and give therefore a straightforward measure of the collective 
hydro velocity at the 'switching' surface T = T c . In Ref. [23] we demonstrated 
that in Pb+Pb collisions at 158 A-GeV the m^-spectra of \24\ can be ex- 
plained simultaneously with the mr-spectra of J/ip and ip' mesons |25| using 
the same set of the hadronization parameters. 

Assuming the fluid freeze-out at constant temperature T, the transverse 
mass spectrum of i-th hadron species in cylindrically symmetric and longitu- 



dinally boost invariant fluid expansion equals approximately to (see Ref. |23 
for further references and details): 



mTdrriT 
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where rriT 



[m: 



+ Pt) 1//2 an d vt is the average transverse flow velocity. 
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Fig. 3: The mr-spectra measured at midrapidity in Pb+Pb at 158 A-GeV by WA97 [g4| for 
f2 and by NA50 |p5| for J/ip and ip' are presented (in arbitrary units) versus tut — m. The 
solid lines correspond to Eq.(|l4|) with T = 170 MeV and vt = 0.19. 

In Fig. 3 we present the fit with Eq.(|14D of the measured my-spectra. The 
temperature is fixed as T = 170 MeV and vt is considered as a free parameter. 
The shapes of all spectra are simultaneously reproduced at vt = 0.19 ± 0.02. 

2 The corresponding calculation for charmonia are not yet performed within this model. 



6 Conclusions 



Statistical hadronization of the QGP is probably an important source of J/ip 
production f5|. Within the SCM || |7| the NA50 data on the J/ip production 
in Pb+Pb at 158 A-GeV can be fitted || (see Fig. 1) for central collisions 
N p > 100. A large enhancement of the open charm over an extrapolation 
from the p+p data is however required 0. A direct measurement of the 
open charm in Pb+Pb at CERN SPS by NA60 would allow to test the SCM 
selfconsistency. 

The SCM predicts that the J/ip suppression at the SPS should be changed 
into the J/ip enhancement in central Au+Au collisions at the RHIC energies 
(see Fig. 2). 

The shapes of m^-spectra for Q, J/ip and ip' are simultaneously reproduced 
(see Fig. 3) in the hydrodynamical picture of the QGP hadronization with 
T = 170 MeV, vt = 0.19 ± 0.02. This supports the hypothesis that formation 
and the kinetic freeze-out of charmonia occurs at the hadronization ]2T)|. 
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